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Abstract 

We investigate the symmetry component of the center variety of polynomial differential 

^\i . systems, corresponding to systems with an axis of symmetry in the real plane. We give 

K^ , a general algorithm to find this irreducible subvariety and compute its dimension. We 

\^ \ show that our methods provide a simple way to compute the radical of the ideal 

^p ' generated by the focus quantities and, therefore, to estimate the cyclicity of a center 

in the case when the ideal is radical. In particular, we use our methods to get a simple 



C^ . proof of the famous Bautin theorem on the cyclicity of the quadratic system. 

o 



1. Introduction 



n ■ We start with a brief history and background. For more information and proofs 

>■ ! for the facts we state without proofs, see either the indicated original source or 



X 



m, 26; 28; m M 



;-( I In his Memoire (|T^, Poincare initiated the study of dynamical systems by 

studying real polynomial differential systems of the form 

I ^ £/(....). (1) 

where U and V are polynomials over the real numbers with a singular point 
(mo, fo), taken to be (0, 0), without loss of generality. There he also defined many 
now standard concepts, in particular, the notion of a center of a system. The 
origin is a center if there exists a neighborhood W of (0, 0) such that every point 
of W other than (0, 0) is nonsingular, and the integral curve passing through that 
point is closed. Moreover, he proved the following theorem. 
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Theorem 1: Assume that the linearization of system (jlp at the origin has 
purely imaginary eigenvalues. Then, without loss of generality, it is of the form 

— = -u + U,iu,v) = U{u,v), (2) 

dv 

— = v + Vi{u,v) = V{u,v), 

where Ui and Vi are polynomials with nonlinear terms. Then the origin is a 
center if and only if there exists a formal power series \1/('U, f ) = -u^ + f ^ + . . ., 
convergent in a neighborhood of the origin, such that 

—u + —v = o. 

ou ov 

The function \l/(M,f) is a local first integral of system (0). 

Lyapunov (|l^ generalized and proved the above theorem for the case when 
P and Q are real analytic functions. 

There are many different ways to enclose the set of plane real system (RS), 
into the set of two dimensional complex systems (CS). The most convenient 
and commonly used way is the following. Consider the real plane {u, v) as the 
complex line with the variable x: 

X = u + iv. (3) 

Then System (|l]) is equivalent to the equation 

^^ = ^(a;,x), (4) 

where P = iU — V. In many cases it is convenient to use just equation (^, 
however it is natural to add to this equation its complex conjugate, —ix = 
P{x, x) and consider y = x a.s a. new variable and Q = P as a. new function. As 
a result, we get the system of two complex differential equations 

ix = P(x, y), -iy = Q{x, y). (5) 

Thus there is one-to-one correspondence between systems from (RS) and the 
subset of (CS) consisting of systems of the form (|) where the second equation 
is the complex conjugate to the first one. 

For polynomial systems of the form (0) the procedure above yields the system 

dx 

^~dt " ^ + Pii^^y) = Pi^^y)^ (6) 

dv 
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where Pi and Qi are complex polynomials with nonlinear terms and Qi = Pi. 
After the change of time idr = dt we can write @ in the form 

dx 

— = x + Pi{x,y) = P{x,y), (7) 

--1- = y + Qiix,y) = Q{x,y). 

In (|TUp, Dulac considered system (|^) where Pi and Qi are arbitrary complex 
polynomials with nonlinear terms. Moreover, he gave the following definition for 
a center at the origin. 

Definition 1: System (0j has a center at the origin if there is an analytic first 
integral of the form 

oo s 
s=3 i=0 

where the Vj^s-j are functions in the coefficients of P and Q. 

If (J^) is the complexification of (0) by means of (^) then this definition is in 
agreement with the definition given by Poincare. 

When P and Q are quadratic polynomials, Dulac gave necessary and sufficient 
conditions on the coefficients of P and Q such that system (|^) has a center at 
the origin. Moreover, he asked if one can find necessary and sufficient conditions 
on the coefficients of P and Q (of any given degree) such that system (|^) has a 
center at the origin. This is the so-called center problem. 

In this paper, we present a partial solution to the center problem for polyno- 
mial systems, using methods from computational algebra. 

Any polynomial system of the form (|^ can be written in the form 

^ = ^_ ^ a,,a;^+V = P(x,y), (9) 



■J^ = y- H b,,x^yP+' = -Q{x,y), 



dy 

= y- 

(p,q)es 
where P{x,y),Q{x,y) G C[x, y] and 

S = {{Pi. q^)\p^ + q^>l,^ = l,...J}C {{-1} U N} X N. 

Throughout this paper, N is the set of nonnegative integers. 

We denote by E{a,b){= C^') the parameter space of (^, and by C[a,6] the 
polynomial ring in the variables Qpg, bqp. 

As we have shown above, in the case when 

y = x, bij = (iji, idr = dt (10) 
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system (|^) is equivalent to the system 

^^ = X - Yl ap,^'^'^', (11) 

(p,g)6S 

which has a center or focus at the origin in the real plane {{u,v) \ x = u + iv}, 
where the system can be also written in the form (0) 

u = —V + Ui{u,v), V = u + Vi{u, v). 

In this case we denote the parameter space by E{a). 

Definition 2: We say that System (^ has a center on the set W C E{a,b), if 
for any point (a, b) & W the corresponding system (^ has a center at the origin. 

It is known (|l^) that one can always find a Lyapunov function \1' of the form 
(|) such that 

— P(x, y) + -g-Qix, y) = gn ■ {xyf + (722 ■ {xyf + (733 ■ {xyf + ■ ■ ■ , (12) 

where the gu are polynomials of C[a, h] called focus quantities. Thus, the maximal 
set V C E{a,b), on which System (^) has a center, is the set where all polyno- 
mials gii,i = 1,2,..., vanish, that is, V is the variety of the ideal generated by 
the focus quantities gu. 

Denote by V(/) the variety of the ideal /. 

Definition 3: The set 

V = y{{gii,g22,---,gih---)) 

is called the center variety of system (|^. 

So, for every point in V the corresponding system has a center at the origin in 
the sense that there is a first integral of the form (|^). However, if (a, 6) G V 
and apq = bqp for all (p, q) G S, then such a point corresponds to a real system 
of the form ([TT|) , which then has a topological center at the origin in the plane 
X = u + iv. (For a geometrical interpretation of the center of the complex system 
(I) see, e.g., (|3|).) 

Therefore, given a system of the form (^), the problem of finding the center 
variety (center problem) of the system arises. 

Among the components of the center variety of polynomial system (§) there 
are at least two components which can be found without computing any focus 
quantities: one component consists of Hamiltonian systems, and the other one, 
which we call the symmetry component, corresponds to systems which have an 
axis of symmetry in the real plane x = u + iv. The component of Hamiltonian 
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systems has a simple geometry because it is equal to the intersection of linear 
subspaces of E{a, b). 

To find the symmetry components one can proceed as follows. With system 
(I) we associate the hnear operator 

QiJ \qij \PiJ \PiJ 

where {pm, qm) £ S. Let M. denote the set of all solutions z/ = (z^i, 1/2, ... , ^21) 
with non-negative components of the equation 

m = ( * ) , (14) 

for all A; G N. Obviously, M. is an Abelian monoid. For A; = 0, finding the 
solutions of equation (JT^) is just the standard integer programming problem, 



see, e.g., ( |34| , Sect. 1.4). Our algorithm is an adaptation of the algorithm to solve 
the standard problem. Let C[A^] denote the subalgebra of C[a,6] generated by 
all monomials of the form 

fl'^i n"'^ ■ ■ ■n"'^ m+im+2 . . . If 21 

"'piqi^P2q2 PiQi qiPi qi-iPi-1 "qip\i 

for aX\ V E M.. In order to simplify notation we will abbreviate such a monomial 

by [u] = [ui, . . . , h'2i]- For u E M., let 

i^ = {1^21,1^21-1, ■■■,1^1) ■ 

Furthermore, let 

IM[z/] = M-[z>], RE[z/] = [iy] + [iy]. 

It is shown in (^; ^^ that the focus quantities of system (^ belong to C[A^] 
and have the form 

9kk= Yl 9{u)IM[ul (15) 

L{u)={k,kY 

with (7(^) G Q, A; = 1, 2, . . .. (Similar properties of the focus quantities were also 
obtained in (|; |5B|).) 
Consider the ideal 

Isy^ = {[y]-[y] \veM)(iC[M\, 

called the Sibirsky ideal of (^. The next statement is an obvious generalization 
of Sibirsky's symmetry criteria (^) for a center for real systems to systems of 
the form (|^). 



Proposition 1: The system ^ has a center on the set V(/ 



symj 
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Definition 4: The set \{Isym) is called the symmetry com,ponent of the center 
variety. 



For points (a, h) G ^^{Isym) for which ttpq = bgp the corresponding systems (|TT|) 
have an axis of symmetry (see Section 3). This is the reason for the name "sym- 
metry component." However we should mention that we do not know a proof that 
V(/sym) is indeed a component of the center variety, that is, it is not a proper 
subvariety of an irreducible subvariety of the center variety. However Proposition 
1 implies that y^{Isym) is a subset of the center variety and we shall show below, 
that Y{Isym) is irreducible. Moreover, for all cubic polynomial systems V^{Isym) 
we have investigated, it is indeed a component. Thus we conjecture that for any 
polynomial system of the form (^, \^{Isym) is a component of the center variety. 
In the present paper we give a simple and effective algorithm for computing 
generators for the Sibirsky ideal, using methods from computational algebra. 
This allows us to find a finite set of defining polynomials for the symmetry 
component of the center variety. Moreover, we prove the following results. 

Theorem 2: The Sibirsky ideal Igym is prime in C[a,6]. 

Theorem 3: The dimension of the symmetry component, 'V{Isym), of the center 
variety is equal to I if all coefficients on the right-hand side of system ^ are 
resonant and I + 1 otherwise. 

Recall that the resonant coefficients are the ones which cannot be canceled by 
transformation of the system to the Poincare-Lyapunov normal form, that is, 
the coefficients of the form a^fe, hkk- 

It should be mentioned that a substantial amount of literature is devoted to 
different particular subfamilies of polynomial systems, mainly to systems of the 
second to fifth degree (the bibliography on quadratic systems alone by J. Reyn 
(pop contains approximately 1500 references). However, the Symmetry Compo- 
nent Algorithm in Section 2, and Theorems |^ and |^ are among the very few 
statements known up to now about the whole class of polynomial systems. 

The center problem is closely connected to the cyclicity problem, which is 
sometimes called the local 16th Hilbert problem (0). 



Definition 5: Let Ua^e be the number of limit cycles of system (1J_) in an e- 
neighborhood of the origin. Then we say that the singular point x = of system 
( |7i| j with given coefficients a* G E{a) has cyclicity k with respect to the space 
E{a), if there exist 5q, eo such that for every < e < eo and Q < 5 < 5q 

max Ua^e = k. 

a£Us(a*) 
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The cychcity of the origin of real quadratic systems (that is, system ( pT]) with 
quadratic non-Hnearities) was first investigated by Bautin (j^), and by Sibirsky 
(pop (see also (pSj)) for the system (^) with homogeneous cubic nonlinearities. 
They proved that the cyclicity of the systems is less or equal two, respectively 
four. (If we take into account perturbations by linear terms, then the cyclicity 
is 3 and 5, respectively). 

In this paper we show that Theorem ^ along with other methods from com- 
putational algebra provides an efficient tool to investigate the cyclicity of poly- 
nomial systems in those cases, where the ideal of focus quantities {gu, g22, ■ ■ ■) 
is radical. Moreover, even when the ideal of focus quantities is not radical, we 
use our algorithm along with some standard methods to compute the cyclicity 
of some cubic systems, see (|1^. 

2. An Algorithm for the Symmetry Component 

In this section we give an algorithm to find a finite set of generators for the 
Sibirsky ideal, hence for the symmetry component of the center variety. It works 
for general systems (^. As a corollary, we obtain a Hilbert basis, that is, a finite 
minimal generating set, for the monoid Ai described in the introduction. 
Let 

^ Cii C21 ... Qi Q2 ... C22 Cu 
Ci2 C22 ... Q2 Qi ... C21 Cii 



A 

be a (2 X 2/)-matrix with entries Qj G Z. For each k eN, let 



V{k) = !^u= iu,,...,U2i) EN'' \A-u' = (^l^Y 

The proof of the following lemma is straightforward. 

Lemma 1: Let M.a = UfcgNV(A;). Then M._a is an Ahelian submonoid ofN'^K 
Also, if 

fi = {fii, . . . , fi2i) e Ma, 

then ft = {^2u ■ ■ ■ , /^i) e Ma- 

Let d = 21 and R = C[xi, . . . ,Xd\ be the polynomial ring in d variables. We 
consider the following binomial ideal in R: 

lA={^^-^^\li^MA). 

We will obtain a Hilbert basis for Ma from a Grobner basis of the ideal Ia- 
If the matrix A arises from the coefficients of system (1) as follows: 

A = 



Pi 


P2 ■ 


■ pi 


qi qi-i . 


.. qi 


qi 


q2 ■ 


■ qi 


Pi Pi-i ■ 


. . Pi 
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then J^ is precisely the Sibirsky ideal of (P). 

We will represent /^ as the kernel of a homomorphism of polynomial rings, so 
we can use a standard algorithm to compute a finite generating set for this ideal. 
As a consequence, we obtain an algorithm to compute a finite set of polynomials 
that define the symmetry component of the center variety of (1). 

Let S = C[t^ , tf , yi, . . . , yi\. Define a ring homomorphism : R — > S by 

™ I , „, .C!_i+i,2 .c;-j+i,i 

for i = 1, . . . , /. 

Theorem 4: ker(0) = /^. In particular, /^ is a prime ideal in R. 

Proof: The second statement follows immediately from the first, since S* is a 
domain. 

To prove that ker(0) = J^, we first show that ker(0) is a binomial ideal. We 
can factor the map as follows. Let 

T = k[ti,t2,si,S2,xi,. . . ,Xd,yi,. . . ,yi\, 

and 

(p' : C[xi,...,Xd] — >T 

be defined like 0, except that if q < 0, then it appears as an exponent of the 
variable Sj instead of tj. Then is equal to the composition of 0' followed by 
the projection 

T^T/(tiSi-l,t2S2-l). 

It is straightforward to see that ker(0) = ker(0'). Let 

J={xi-y/;H';\xi+,-yi^i+^fl-'*''H'^-'+''' | z = 1, ...,/) C T. 

Then it follows immediately from (|I|, Theorem 2.4.2) that ker(0) = ker(0') = 
Jri-R. We obtain a generating set for Jfl-R by computing a reduced Grobner basis 
for J using an elimination ordering with Xi < yj, tk, Sr for all i,j, k, r, and then 
intersecting it with R. Since J is generated by binomials, any reduced Grobner 
basis also consists of binomials. This shows that ker(0) is a binomial ideal. 

Now, let x" — x^ be a binomial in R. We may assume that the two monomials 
have no common factors, that is, supp(a;) fl supp(/5) = 0. Then 

0(x" - x^) = {yiti^nrr ■ ■ ■ {yitrtrriyitrtrr^' ■ ■ ■ {yitrtrr 

-ivitT'trf ■ ■ ■ {vitTtTtivitTtTt^' ■ ■ ■ {yitrtrf' 

_ (01+02;) „ ("i+Oi+i) .(ciiOiH hcna;+C;2"i + lH hci2ad) .{ci2QlH hc;2Q!+Ciia;+lH hcii«d) 

— iJl ' ' ' iJl ^1 ^2 

, (P1+P21) ^ (A+A+i) .(cii/3i+-+qi/3i+q2A+i+-+ci2/3d) (ci2/3i+-+q2A+caA+i+---+cii/3d) 
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Thus, 0(x" — x^) = if and only if 

for alH > 1, a, + ad-i+i = f3i + Pd-i+i, 

ciiai H h cnai + Ci2ai+i + h Ci2ad 

= CuPi H h QiA + Q2A+1 H \- CuPd, 

and 

Ci2ai H h Q2«i + Qi«i+i + h ciia^i 

= C12P1 H h C12P1 + cnPi+i H h cii/3d. 

Since supp(a;)nsupp(/3) = 0, we obtain the following facts from the first condition 
in (|T6|). If, for some i, ttj = = /3j, then ad-i+i = = Pd-i+i, so that Pi = ad-i+i- 
Now suppose that /3j 7^ for some i. Then Oj = 0. But then Pd-i+i = 0, otherwise 
ad-i+i = 0, which cannot be since (3i > for all i. And this, in turn, implies that 
Pi = ad-i+i- Hence 

P = {Pi, ...,Pd) = {ad, ...,ai) = a. 

It follows from the last two equations in ([T6|) that a G M.a- This completes the 
proof of the theorem. □ 

Theorem 5: With notation as above, let G be a reduced Grobner basis of I_a, 
with respect to some term ordering on R. Then 

1. The set 

T-L = {^i,fl\yi^ -^ eG}U{e, + ej\i = d-i + l;i = l,...,l}, 

where e^ = (0, . . . , 0, 1, 0, . . . , 0) is the ith basis vector, is a Hilbert basis of 
M.A- Note that e^ + e^-j+i G M.a for all i. 

2. i^= (x^^-x^ IfieH). 



Proof: The proof is similar to that of (p^ Algorithm 1.4.5). It was shown in the 
proof of the previous theorem that any binomial in I^ is of the form x'^ — x^ for 
some /x G Ai^- Hence the Grobner basis G of I^ is of the form 

G = {x^i -x'^, ...,x'^'--x^}. 

We first show that the set 7i is a generating set for A^^. Suppose not, so there 
exists fi G M-A which is not an N-linear combination of elements in 7i. We can 
choose /i so that x'^ is minimal with respect to the chosen term ordering of R. 
Since /x G ?^, we have 0(x'^ — x'^) = 0, so that x^ — x^ G Ia- Hence the leading 
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10 



term of this binomial is divisible by a binomial in G. That is, there exist i, j3 
such that 



x^-x'^ 



x^fx^' 



x^O +x'"'(x' 



^J■^l^P _ Y^ 



x''^ 



Thus x^^fx^^ 



x*^ — x*- 



G I A- But /_4 is a prime ideal, and it is immediate from the 
definition of that it contains no monomials, so (x'^ — x'^) G /^t- Moreover, 
x/3 < x'^. Hence, by the choice of /x, /5 is a linear combination of elements of 7i, 
which implies that /x = /5 + /ij is also a linear combination of elements of 7^. This 
is a contradiction to our assumption on /z. Thus, ?i is a generating set for M.ji,. 
To show that Ti is minimal, suppose that some /i, or fTi is an N-linear com- 
bination of elements in Ti. Since the Grobner basis G is reduced, the linear 
combination cannot contain any summands coming from the fii,JIi. But observe 
that all the vectors e^ + e^ are symmetric, so that 



X 



X 



0. 



A similar argument disposes of the other cases. Thus, we have shown that H is 
the Hilbert basis of Ma- 
2. {^>'-xP\nen) = {G) = Ia- □ 

We can also compute the dimension of the afiine variety V(/^) of the ideal 
Ia- 

Theorem 6: The dimension of\{lA) is equal to I if Cn = Q2 for alii = !,...,!, 
and I + 1 otherwise. 



Proof: According to (|33| , Lemma 4.2) the dimension of V(/^) is equal to the 
number of linearly independent column vectors in the matrix 



B 



( 1 

1 



Cll C21 

\ C12 C22 







1 1 

Qi Q2 

Q2 Qi 



1 \ 

1 



C22 C12 

C2I Cll / 



/I 


... 





1 ... 





... 





... 


V 


... 













1 1 

Q2 - Qi 






1 



1 







C22 ~ C21 C12 — Cll 

/ 



The theorem now follows. 



D 
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Proof of Theorems ^ and |^. Theorem |^ is a corohary of Theorems ^ and ^ and 
Theorem |^ follows immediately from Theorem p. 

To close this section, we summarize the algorithm, and in the next section we 
shall provide several examples. 

Symmetry Component Algorithm 

Input: Two sequences of integers Pi, . . . ,Pi,Pi > —1; qi, . . . , qi, qi > 0. (These 
are the coefficient labels for a system of the form @).) 

Output: A finite set of generators for the Sibirsky ideal Igym of (^, and the 
Hilbert basis Ti of the monoid M.. 

1. Compute a reduced Grobner basis G for the ideal 

in C[a, 6, 2/1, ... , yi, tf , t^], with respect to any elimination ordering with 
{ti,t2} > {yi,...,yd] > {apig,,...,6g,pj. 

2. Isym={GnC[a,b]). 

3. H = {u, u \ [u] — [u] e 0} Li {ei + Bj \ j = d — i + I; i = l^ . . . , 1} is a Hilbert 
basis for the monoid A^, where ej = (0, . . . , 0, 1, 0, . . . , 0) is the ith basis 
vector. 

3. Invariants of the Rotation Group and the Symmetry 
Components 

The symmetry component of the center variety of a real system in the com- 
plex form (|TTD was investigated by Sibirsky and his coworkers (|3T| ; |32D , using 
invariants of the rotation group of this system. 
He showed that a monomial 

piqi P2q2 PlQl QlPl q2P2 qipi 

is invariant under the action of the rotation group 

X I — > xe''^ (16) 

if and only if z/ G N^' is a solution of the Diophantine equation 

(Pi - qi)i^i + {P2 - q2)i'2 + ^{Pi- qi)^i (17) 

+ {Qi - Pi)t^i+i ^ h (gi -pi)z/2« = 0, 



which is obtained by subtracting the second equation of system ( [14D from the 
first one. 

Proposition 2: The set of all non-negative integer solutions of equation (|7^ 
coincides with the monoid M. defined by equation j^) . 
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Proof: Obviously, every solution of (|1^ is also a solution of (jl^). Conversely, let 
z/ be a solution of ([17|) . Then 

L\u)=L\u) = k (18) 

and 

L\e,) + L'iei) = L\e2i-i) + L^e^i^i) =Pi + qi>l, (19) 

for i = 1, . . . ,1. Taking into account that -L(z/) is a linear operator, we conclude 
that k on the right-hand side of (|18D is non-negative. □ 

Thus, to find invariants of system ( pA]) under the rotation group (|T6|) it is 
sufficient to find a generating set of the monoid A4. It also is easily seen that 
the monomials [i'],^' G Ai are invariants of the system (|^) under the action of 
the transformation 

X I — > xe*"^, y I — > ye~^'^. 



Knowledge of the invariants of system (|TT]) allows one to determine whether 
the corresponding vector field has an axis of symmetry. Denote by E the subset 
of E{a) consisting of systems ([TT|) such that the fraction U{u, v)/V{u, v) is irre- 
ducible, where U, V are the functions on the right-hand side of system (0). The 



following theorems are proven in (|3lD . 



Theorem 7: If a* e E{a) and IM[i/] = Q at the point (a*, a*) for all v E M, 
then the corresponding vector field (j^ has an axis of symmetry passing through 
the origin. Moreover, if a* E E then the opposite statement holds. 

According to this theorem, for every point in V(/sym) of the form {a*, a*) the 
corresponding vector field (0) has an axis of symmetry passing through the origin. 
This is the reason for calling \^{Isym) the "symmetry component." However, we 
have not investigated whether for arbitrary points in V^{Isym) the corresponding 
systems (^ have any kind of symmetry. 

As an immediate corollary of Theorem ^ we get 

Theorem 8: If a* e E{a) and IM [v] = at the point {a*, a*) for all u E M, 
then the corresponding vector field (^ has a center at the origin. 

Similarly, it follows from Proposition ^ that if 

IM[z/] = (20) 

for all u E M at {a*,b*) E E{a,b), then the corresponding system (^ has a 
center at the origin. 

As mentioned above, according to ( PT| ; ^2] ) the monomial [u] is an invariant of 



the rotation group ([T6| ) if and only if t/ is a solution of (p!?!). Proposition 2 implies 
that the set of solutions of equation ( p!7D coincides with Ai. The invariant [z/] is 
called irreducible if u cannot be written in the form u = fi + 6, where /i, 6 are 
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invariants. Therefore the set {[I'i]} is the set of all irreducible invariants if and 
only if {z/j} is the Hilbert basis of A^. 

It is proven in (|3^), that the maximal degree of irreducible invariants of the 
group (1^) is less than or equal to 

2(1+ max (p + q)). (21) 

{p,q)eS 

Using this and equation (|T3|) or (^7\) it is possible to find a generating set of the 
monoid Ai by sorting. In particular, in (^) center conditions are obtained for 
the cubic system 

i— = x(l — aiox — aoix — a_i2X~ z —0202: — auxx — aQ2X — a_i3a;"x) (22) 

UjL 

which can also be obtained from Theorem |^ below by taking imaginary parts. 

Theorem |^ provides another way to find a Hilbert basis of the monoid A4 
and, therefore, the center conditions (^). We apply it to find the symmetry 
component of the general cubic system 

X = x{l- aiox - aQiy - a_i2X^^y'^ - a2ox'^ - auxy - ao2y^ - a_i3x"-^y^), (23) 
y = -y{l - b2-ix^y'^ - hox - boiy - bs-ix^y'^ - b2ox^ - buxy - &02l/^), 

where x,y,aij, b^ G C. 

Theorem 9: The symmetry component Vgym of the center variety of cubic sys- 
tem (^^ is defined by the following equations: 

= Oil — bii = aoi^02^2,-l ~ 6^-1,2^10^^20 = ^01002^2,-1 "~ O_i^2^10^20 

= a]^Qa_i3 — bs^-ibfji = aioa-12^20 ^ &oi&2,-icto2 = flioct-i2^io ^ ctoi^2,-i&oi 

= 020^02 ~ ^20^02 = Ctlo'^-l,2&10 ~ ^^01^2,-1^01 ~ '^10^02&10 ~ ^^01^^20^01 

= 00162,-1 — d-ubiQ = aioO'02&io ~ ^01^20^01 = cii^a-u — &2,-i&oi 

= aioa_l,3&2,-l ~ Ct-12&3,-1^01 = '220'J^-1,3^20 ~ Ct02&3,-1^02 = '2io&02 ~ '^20^01 
= '^'-02^3,-1 — 0,-13^20 — Q'0lO-12&3,-l ~ Ct-13^2,-l&10 = Ooife20 ~ O02&10 
= a2oa_i3 — ^3,-1^02 ~ C^lO'^-13^20^10 ~ Ooiao2&3,-1^01 
= aioa20'^-l,3&10 — Ct01&3, -1^02^01 = ^^10^02^2,-1 — Ct-12a.20&01 

= aiQao2 — ^20^01 ~ Cfl0O02&02&2,-l — Ct-12a20^20^01 (24) 

'02&10 



'^10^^02^2,-1 ~ 0-12^20^01 — ^^01^3,-1^02 — 020^^-13^10 ~ ^0-^020 — &02& 



0010-12^20 ~ Oo2&2,-1^10 = OiqO_i3620 ~ Oo2&3,-1^01 ~ O_i2O20&10 — ^02^2,-1^01 
O0lO-12O20^20 ~ Oo2&02&2,-l&10 = OqiOo2&3,-1 ~ O_i3620^10 

O10O20O-I3 ~ ^3,-1^02^01 ~ OoiO„i202o " ^02^2,-1^10 = OioO_i36io ~ OQ163 _i6oi 
O10O-I3&IO ~ '^01^3,-1^01 ~ OiqO_i36io — 00163,-1^01 = Ooi&3,-1 " 0_i36iq 
OoiO-l3&20^2,-l — O_i20o2&3,-1^10 = 0010200-1362,-1 — 0-1263,-1602610 
0100-1263,-1602 — 0200-1362 1601 = O10O-12O0263 1 — 0-1362062,-1601 
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'^-12^3,-1^20 ~ fl02O-13&2 -1 = Cf-12'^20^3,-1 ~ Ct-13^02&2 -1 



Q-loO-12^3, 


-1^10 ~ CtOlC^-13^2,-1^01 = CtoiCt-13^2,-1 ~ '^-12^3,-1^10 


'^-12^20 ~ 


'^02^2,-1 = '3'-12'^20&20 ~ '^02^02^2,-1 = '^-12'^20^20 — '^02^02^2,-1 


«?0«-12^3, 


-1 — '^-13^2,-1^01 = '^-12'^20 ~ ^02^2,-1 = '^-12^3,-1^02 " '2200-13^2, 


«'-12&3,-l 


— a„i3fe2,-l = '3--i2'^02^3,_i ~ '^-13^20^2,-1 ~ '^lO'J^Ol ~ ^10^01 


2 l3 

o.oia._i302 


_1 — a_i2&3,-1^10 = Ctioa_i2&3 -1 — O-_i3&2,-1^01- 



Proof: It is enough to show that the above equations (^J) form a Grohner basis 
of the ideal Isym- To compute Isym, one can use the Symmetry Component Al- 
gorithm, with any computer algebra system. Here we are using the specialized 
system Macaulay ( p!3| ) which performs Grobner basis calculations substantially 
faster than most general purpose symbolic calculation packages. Figure |l] shows 
the Macaulay session used to compute Igym for system (p3|). To simplify nota- 
tion, we renamed the variables aij and bij as follows: Xi = aio,a^2 = '^^0152:3 = 
a_i2, X4 = a2o, x^ = ail, xq = ao2, x-j = a-1,3, xg = &3,-i; xg = 620? 2:10 = &11, a;ii = 
bo2, a;i2 = &2,-i, a;i3 = &10, xu = &oi- n 

Recently some sufficient center conditions were obtained for the real systems 
with homogeneous nonlinearities of fourth and fifth degrees (^; |^). The following 
are the general complex forms of such systems. 



X = x{l- asox^ - a2ix^y - ai2xy'^ - a^-^y^ - a^^x ^y^), (25) 

y = -y{l - bi-ix'^y'^ - b^ox^ - b2ix^y - b^xy^ - 6o,32/^), 



and 



X = x(l — a4ox'^ — a^ix^y — a22x'^y'^ — ai3xy^ — ao^y^ — a-i^x ^y^), (26) 
y = -y{l - b^^ix^y'^ - b^^x^ - b^ix^y - b22X^y'^ - bi^xy'^ - b^^y'^). 

The following theorems give some center conditions for complex systems with 
homogeneous nonlinearities of fourth and fifth degrees, respectively. 

Theorem 10: The symmetry component of the center variety of system i\2di ) is 
defined by the following equations 

= 030(3.03 — ^30^03 = 021012 — &21&12 = 030^12 ~ '^21^03 

= a3oa.l2&i2 ~ '^21^21^03 = '330'J^12^12 ~ '321^21^03 = '32l'303^21 ~ '^12 ~ ^30^12 

= '32l'203^21 — '312^30^12 = '^^12^30 ~ '^03^21 = '22l'^03 ~ ^30^12 

= 030^12 — ^21^03 = O30'3-l,4&30&12 " O2lO03^4,-l&03 

= a3Qa_i^46i2 — 0.21^4,-1^03 ~ 'J^30'^-l,4&21 ~ '^12^4,-1^03 

= a2ia-l,4&30&21 — ai2ao3&4,-l&12 = O30a2lO_l,4&21 — a.i264,-l&12&03 
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il : K = ZZ/31991; 

12 : R = K [y_l . .y_7, s, t, n, m, MonomialOrder -> Lex] ; 

13 : I = Ideal (s*n-l, t*m-l); 
o3 : Ideal of R 

14 : Rbar = R/ {I*R) ; 

15 : S = K[c,x_l . .x_14, MonomialOrder => Lex]; 

16 : f - map (Rbar, S, matrix {Rbar, { {n*in, y_l*s, y_2*t, y_3*n*t^2, y_4*s'^2, y_5*3*t, y_6*t''2, 

y_7*n*t''3, y_7*m*s^3, y_6*s'^2, y_5*s*t, y_4*t^2, y_3*m*s'^2, y_2*s, y_l*t } } ) ) ; 

06 : RlngMap Rbar < S 

n : J = ker f 

07 = ideal (x -x ,xx -x x ,xx -xx ,xx x -xxx ,xxx -xxx , 

5 10 12 13 14 4 6 9 11 1 11 13 2 4 14 1 6 13 2 9 14 

2 22 22 22 2 

XX -XX XX -XX ,xx -XX ,xx -x x ,xx X -xxx , 

1 11 4 14 2 9 6 13 16 9 14 2 4 11 13 2 11 12 3 4 13 

2 2 
xxx -xxx ,xxx -XX X ,X X -X X X ,xxx -XX X , 

2 6 12 3 9 13 13 9 6 12 14 13 4 11 12 14 1 3 13 2 12 14 

2 2 3 3 3 3 

XXX -XX X ,xx -XX ,xx -X X ,xxx -xxx ,XXX -XXX , 

1 3 13 2 12 14 2 12 3 13 13 12 14 1 7 12 3 8 14 4 7 9 6 8 11 

2 2 2 2 

XX -XX, XXX -XX X ,xx -XX ,xxxx -xxxx ,XXXX -XXX X , 

68 79 238 7 12 13 47 8 11 179 13 268 14 147 13 23 11 14 

2 2 2 2 2 2 

XX X -xxx ,xxx X -xxxx ,XXX -XXX ,XXX -XXX , 

1 11 12 3 4 14 1 6 11 12 3 4 9 14 1 6 12 3 9 14 2 8 11 4 7 13 

2 2 2 2 2 2 

XXX -XX x,xxx -xxx ,XXXX -XX X X ,XXX -XXX , 

239 6 12 3 179 68 14 2349 6 11 12 13 268 79 13 

2 2 22 33 2222 

XXX -XX X ,xxx -X X X ,XX X -XXX ,XXX -XXX , 

14 7 8 11 14 2 3 4 11 12 13 1 7 13 2 8 14 1 7 13 2 B 14 

3 3 4 4 4 4 

XXX -xxx ,XX -XX ,XX -XX ,XXXX -XXXX ,XXXX -XXX X , 

17 13 28 14 28 7 13 17 8 14 279 12 368 13 247 12 38 11 13 

2 2 2 2 

XXXX -XXX X ,xxxx -xxx X ,XXX -XXX ,XXX -XX X , 

138 11 47 12 14 1368 79 12 14 389 67 12 348 7 11 12 

2 2 22 2223322222 

XXXX -XXX X ,XXX -XXX ,XX -XX ,XXX -XX X , 

1 3 8 13 2 7 12 14 2 7 12 3 8 13 3 9 6 12 3 4 9 6 11 12 

22 22 22 22 23 32 22 22 

XXX -XX X ,xxx-xx X ,xx -X X ,xxx -xxx , 

3 4 9 6 11 12 13 8 7 12 14 3 4 11 12 3 8 11 4 7 12 

22 22 23 32 32 23 43 34 

xxx -xxx ,xxx -xxx ,XXX -XX X ,xx -XX ) 

368 79 12 27 12 38 13 138 7 12 14 38 7 12 

o7 : Ideal of S 



Figure 1: Macaulay computation for system(23) 



Ctl2Ct-l,4&30 ^ '^03^4,-l&21 = a.3oa.l2Ct-l,4&30 — Ct03&4,-l&21&03 

0210-1,4^30 ~ ^^03^4,-1^1, 2 — 'J^2lOo3&4,-l "~ 0'-l,Ab3Qbi2 

O'30O'l2O'-l,4: — ^4,-l&21&03 ~ O3oa2l'^-1,5 — ^4,-1^12^03 (27) 

Ct2l'2— 1 4^21 — (l^2^A: —1^12 ^ Ct2l'-^— 1 4^21 — (1^2^^ —1^12 

0210-1,4^21 ~ '2l2^4,-l&l2 = 0210-1,4^21 — Oi2&4,-l&l2 

'^12^4,-1 ~ 0_i,462i = 02iO_i,4 — &4,-l&12 

O30O_i 4630621 — O12OQ364 __i6o3 = 03qO_i 4630621 ~ O12O0364 _i6q3 

%0'^-l,4^21 ~ 01264^603 = 02lO_i 463Q — OQ364 1612 
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O30Cf2lO_i,4&30 ~ Oo3&4__i&12&03 — 03o'^2lO_i,4&30 ~ '203&4,-l ^12^03 

ai2ao3^4,-l ~ '^_l,4&30^21 = Ct3o'^2lCt-l,4 — ^4,-1^12^03 

^300-1,4030 ~ O03''4,-l003 — ^300-1,4030 - a-0304,-lOo3 
2 !,3 1,3 4 ^3 t ^ t3 xA 






30"'-l,4"30 ~ O03&4 1&03 — O3oa_i 4^30 — 0-03^4,-1^03 



flS J,3 _ 3 
"03"4,-l "- 



1,4 



O30 



-6 



'^30*^-1,4 O4 i_6q3 



Proof: Similar to the proof of Theorem ^ To simphfy notation we rename the 
variables aij and hif xi = 030,2:2 = 02152:3 = 012,3^4 = Oo3)2;5 = a-1,4, Xe = 
^4,-1)2:7 = &3o,a;8 = b2i,XQ = 612,3^10 = ^03- Figure | is the Macaulay session 
used to find L,,m- n 



sym- 



12 : R ^ K[y_l . .y_5, s, t,m,n, MonomialOrder ^> Lex] ; 

13 : I ^ Ideal (s*m-l,t*n-l} ; 
o3 : Ideal of R 

14 : Rbar = R/(I*R); 

15 : S = K[c,x_l. .x_10, MonomialOrder => Lex] ,• 
16 : f ^ map(Rbar, S, matrix{Rbar, { {n*m, y_l*s^3,y_2*s^2*t,y_3*s*t^2, y_4*t^3, y_5*m*t^4, 

y_5*s'4*n,y_4*s''3,y_3*s"2*t,y_2*3*t"2,y_l*t'3}})); 

06 : RingMap Rbar < — S 
n : J = ker f 

3 3 2 2 2 2 2 

07 ^ ideal (xx -xx ,xx -xx,xx -xx ,xxx -xxx ,xxx -xxx ,xxx 

14 7 10 2 3 8 9 19 2 10 13 9 2 8 10 13 9 2 8 10 2 4 8 
o7 : ideal of S 



Figure 2: Macaulay computation for system(25) 



Theorem 1 1 : The symmetry component of the center variety of system i\^) is 
defined by the following equations 



= a22 — O22 = 040004 ■" O40604 = 03iai3 — 631 Oi3 = O4o6;^3 — O3j^0o4 

O13&40 ~ O04&31 

O31O04 ~" ^40^13 ~ O4oa;^3 — 631604 = 0490-1,5631 — 01365,-1604 



O22 ~ 622 — O4oao4 — 640604 — 03iai3 — 631613 

O40O13613 — O31631604 = 031004631 — 0.13640613 

,2 „ T ,2 2 t2 



0310-1,5640 ~ 00465,-1613 = 01300465,-1 — 0-1,5640631 

O40O3lO-l,5 — 65,-1613604 = 0400-1,5640013 — 03100465,-1604 
0400-1,5613 — 03165,-1604 = 0310-1,563]^ — 0^^365,-16 



13 



^31 



0-1,5631 ~" Oi365,-i6;^3 — 0130-1,5640 — 00465,-16 



31 



0-1,5649613 



0400130-1,5640 — 00465,-1631604 — 031004265,-1 

'^-1,5631 = 04qOi30_i,5 — 65,-1631604 = 03-^0-1,5 



(28) 



'^13^5,-1 

2 J 2 2 J 2 

O40O-I 5O40 - O04O5 



65,-16 



3 
13 



-ib~ 



04 



^40 



0_;^ 5640 — O0465 -l6| 



2 
04 



n-^ 6^ 

O04"5, 



fl^ 6^^ 

"-1,5"20 



'^40*^-1,5 



"l 



J>. 



04- 
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Proof: Similar to the proof of Theorem ^ □ 

Thus we have presented an efficient algorithm to compute the symmetry com- 
ponent of the center variety. Up to now, the only known method for finding this 
component is due to Sibirsky. His algorithm is as follows. 

1. He gives the formula (^) for an upper bound for the degrees of the irre- 
ducible invariants. 

2. With this bound one can find all irreducible invariants by sorting. 

In ( |36D a method is given to find all "elementary Lie invariants" (our Hilbert 
basis). There the problem is reduced to finding non- negative solutions of a Dio- 
phantine equation similar to our equation (|l^, which is done on a case- by-case 
basis by inspection. 

Moreover neither Sibirsky nor Yi-Rong & Ji-Bin have an analog of our The- 
orem 1, which shows that the obtained invariants generate a prime ideal. But, 
as we will see in the next section, this fact is an important characterization of 
this subvariety of the center variety and is very helpful in the investigation of 
the cyclicity problem. 

4. Applications to Cubic Systems 

In recent years many studies have been devoted to investigating different sub- 



families of the cubic system (|23| ) (see, e.g., (§; |; ||; |ri|; ^ |3g; ^) and references 
therein). In the case when conditions ( [T0| ) are satisfied system (|23| ) is equivalent 
to system (^2]). It should be mentioned that the center- focus problem is much 



better investigated for the real cubic system (E^T) than for the general system 



In this section we will show that the results obtained above, together with ad- 
ditional tools from computational algebra, give a simple efficient way to compute 
the radical of the ideal of focus quantities and, therefore, to solve the cyclicity 
problem in those cases where this ideal is radical. As a rule it is easy to find 
rational parameterizations of components of center varieties. This fact, together 
with the following theorem and Theorem ^ gives an easy method of finding the 
radical of the ideal of focus quantities. For general rational parameterizations 
this is a difficult computational problem. 

Theorem 12: // the variety V( J) of an ideal J of C[xi, ... , Xn] admits a ratio- 
nal parameterization 

Xi = — r-, z = l,...,n, and 

diyti, . . . ,tm) 

C[xi, ...,Xn]n{l- tg, gi{ti, ..., tm)xi - fiiti, . . . ,tm) ■■ i = I, ■ ■ ■ ,n) = J 
(where g = gig2 ■ ■ ■ gn), then the ideal J is a prime ideal o/ C[a;i, . . . , x„] . 
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Proof: It is sufficient to show that the ideal 

H = {1 - tg, giiti, ..., tm)xi - fi(ti,...,tm) ■■ i = l,---,n) 
is prime in C[xi, . . . , Xn, ti, . . . , tm, t]. Consider the ring homomorphism 

^ : C[a;i, . . . , x„, ti, . . . , tm, t] — > C(ti, . . . , t^), 
defined by 

ii i 

9i 9 

i = 1, . . . ,n,l = 1, . . . ,m. It is sufficient to prove that ker(-?/') = H. It is clear 
that H C ker (■?/;). We will show the other inclusion by induction. 
Let us suppose that h G ker^ip) and h is linear in x, that is, 

n 

h = ^ «j(ti, ...,tm, t)xi + ao(ti, . . . , tm). 
Then 

n 

c^otg = -"^^ifitgi, 

i=l 

where gi = g/ gi- Therefore 

n 

h = ^ ai{xigi - fi)tgi + {l-tg)h. 

i=l 

Hence, h E H. 

Assume now that for all polynomials of degree k in xi, . . . ,Xn and h G ker (■?/'), 
we have h & H. 

Let h G ker(-?/') be of degree k + 1 in xi, . . . , a;„. We can write h in the form 

n 
n / ^ i^iyXi, Xj-)_i, . . . , Xn, ti, . . . , Zm, Z) ~r iloyZl, ■ ■ ■ , Zm, ZJ 
i=l 

(here every term of hi contains Xj). Consider the polynomial 

u = Y] —{xigi - fi)tgi + (1 - tg)h. 
1=1 

Then u = h + v^ where 



h 

V = -ty^^figi— -tgho. 

■ , -^J 

J=l 

Since h,u & ker(-?/') we get that v G ker(-?/^). Then, by the induction hypothesis, 
V & H and, hence, h & H . □ 
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We now apply the results obtained so far to the investigation of the cyclicity 
problems in some specific cases. First consider the systems with homogeneous 
quadratic and cubic nonlinearities, 

X = x- aiox'^ - aoixy - a^^y'^, ,^^. 

y = -{y -bioxy -boiy^ -b2,-ix^) 

and 

X = X- a2ox^ - aiix'^y - ao2xy'^ - a_i3?/^, /^qx 

y = -{y -bo2y'^ -bnxy^ -b2ox^y -bs-ix^), 

correspondingly. 

In the case when the conditions (|TD|) hold and the linear perturbations also 
are taking into account the system p9| ) corresponds to the real system on the 
plane {u,v), x = u + iv 

ix = iXx — X — aiox^ — aoixx — a-ux"^ (31) 

and the system (|30|) corresponds to 

ix = iXx — X — a2ox'^ — anx^x — ao2Xx'^ — a-iax^, (32) 

where A G M. 

For the first time the cyclicity of the origin of the system (|3T|) was investigated 
by Bautin (||) and of the system ( ^2]) by Sibirsky (^) (later on another proofs 
were obtained by Zol§dek (|37| ; [38|) and Yakovenko (|35D). They proved that the 
following statement holds. 

Theorem 13: The cyclicity of the origin of system ^ ) equals 3 (0; ^; |||j 
and the cyclicity of the origin of system ^3^ ) equals 5 ((5^; |M 



The crucial and the most difficult part of the proofs of this theorem is the 
following statement (see e.g. (|^; ^ ^ |3^ for detail derivation Theorem ^ 
from Proposition ^). 

Proposition 3: 1) The first three focus quantities generate the ideal of focus 
quantities of systems (^^ and P% ) with A = 0. 

2) The first five focus quantities generate the ideal of focus quantities of sys- 
tems ii5D) and iWE) with A = 0. 



Bautin proved Proposition ^ for the real quadratic system (pTD in the Kapteyn 
form. The proof is quite complicated, because the ideal of focus quantities is not 
radical in this case. A simpler proof on Bautin's way was given by Yakovenko 

Zol§dek (0; ^S|) found a new way. He proved the Proposition ^ for the ring 
of polynomials, which are invariant under the action of the rotation group and 
showed that the ideal of focus quantities is radical in this ring. 

We will show that with Theorems 1 and O the proof of the celebrated Bautin 
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theorem, as well as the treatment of the cyclicity problem for the system (pOD, 
becomes straightforward, using only basic knowledge of computer algebra. 

Consider first the system (^). Let I^'^^ = {gii,g22, . . .) be the ideal generated 
by all focus quantities of this system (the so-called Bautin ideal) and let l\^ — 
{gii-i 922, ■ ■ ■ , Qkk) be the ideal generated by the first k focus quantities. 

Computing the first five focus quantities by means of the algorithm given in 
(pl|), and then reducing them, we find that 

fl'ii = c^ii — ^ii; 

5'22 = O'2qO'02 — &02&20; 

5'33 = (3a20'^-13 + 8O20'^-13^20 + 3ao2^3,-l 

- 8ao2&o2&3,-i - 3a_i3&20 " ^^mh-i)/^] 
g44 = (-9a2oa-i3&ii + aua^isbl^^ + 9aii6o2&3 -1 - cto2^ii^3 -i)/16; 
955 = (-9a2oa-i3&02&2o + a2oao2a-i3&2o + 9a2oao2&o2^3 -1 + 

18a2oa^i3&2o&3 -1 + 60020-13^3 -1 - '^02 + 

&02&20&3,-l - 18ao2a-13&02&3 -1 - 6a^l3^20^3 -i)/36. 

Furthermore, the variety of the ideal / coincides with the variety of the first five 
focus quantities. More precisely, the following result holds. 



Theorem 14: The center variety \ {I ^'^^) of the system (UW consists of the three 
irreducible components: 

V(/(^)) = V(jf ) = V(JO U V(J2) U V(J3), (33) 

where Ji = (an - &ii,3a2o - &20,3&02 - ^02), J2 = (aii,&ii,a2o + 3620,^02 + 
3ao2, a_i363^_i— 4ao2&2o) and J^ = (0200-13 ~ ^3,-1^02' ^20^02 ~&20^02; 0200-13620" 

002^3,-1^025 On — 611, Oo2&3,-1 ~ O_i3620/- 



Proof: It is easy to check (using, for example, the radical membership test, see 
e.g. (D) that V(4'^) = V( Ji) U V( J2) U V( J3). Thus we only have to show that 
5'mm|v(jfe) = for all m > and k = 1,2, 3. 

Indeed, the systems corresponding to the points of V( Ji) are Hamiltonian. 

For the variety V(J2) in the case a_i3 = (4ao2&2o)/63,-i one can easily find an 
invariant conic h and invariant cubic I2 and check that the system has the first 
integral 

$ = Iflf, (34) 

defined on V(J2)\V(i7), where H = (&20&3 -1, ^l -i(oo2&3 _i + 462o))- Computing 
the quotient of the ideals (e.g. by means of the algorithm from (^)) we get 
J2 : H = J2. Therefore, 

V(J2)\V(/J)=V(J2). 

This implies that the system (|30D has a center on the whole component V(J2). 
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Finally, according to Theorem ^ ^{Js) is the symmetry component of the 
center variety. 

The irreducibility of the components V( Ji) and V( J2) is obvious, and V( J3) 
is irreducible because due to Theorems |^ and ^ the ideal J3 is prime. □ 



A very similar theorem is proven also in ( pq ) (and in (p!2| |38D for the real case). 
However, our theorem also shows that ( |33D is the irreducible decomposition of 
the center variety. 

It is worthwhile to mention that, if one presents the solution of the center 
problem for a system of the type (|^) writing out the irreducible decomposition 
of a center variety, then the answer is unique. Otherwise it can happen that the 
center conditions for the same system obtained by different authors can look 
very different. In fact the center conditions for system (|30|) were first obtained 
by Sadovsky (pT]). He found eight such sets but they look very different from 



the ones presented in Theorem |T4[ However computing the Zariski closure of the 
sets given by Sadovsky and taking intersections of the corresponding ideals using 
standard algorithms from computational algebra (see below) we get the ideal J5 . 
This implies that his center conditions coincide with the three components given 



in Theorem 14 



Recall that we can use the following algorithm to compute the intersection of 
two ideals (|^). Let / = (/i, . . . , /,.) and J = {hi, . . . ,hs) he ideals in k[xi, . . . , x„]. 
Compute a Grobner basis for the ideal 

(t/l, . . . , tfr, (1 - t)hi, . . . , (1 - t)hs) C k[t, Xi,...,Xn] 

using a lexicographic term order with t greater than the Xi. Those elements of 
this basis which do not contain the variable t will form a basis of / fl J. 

For the case of quadratic system one can easily compute the three first focus 
quantities and get 

5^22 — '^10^^-12^10 "'^01^01^2, -1^3(^^-12^10 ^'^01^2, -1)^3(^01^01^2, -l—Ctlo'^-12^1o)) 
5^33 = 12(^01*^10^01^2,-1 — fiio'^-12&01^lo) ^ 48('2oi^01^10^2,-l — CtOl'J^10'^-12&io) 

+ 8(^^01^^-12^01^2,-1 ~ '^10'^-12^10^2,-l) ^ ig ('2oi'^-12&01^2,-l ~ '^10'^-12^10^2,-l)- 

Here and below the upper index (g) means that we are speaking about the 
quadratic system. 



Using the equations above and Theorems 1 and |T2| it is easy to see (similarly to 
the case of the system (^)) that the following theorem holds (see ( p!0| ; |12|; ^; ^7^ 
for details). 



Theorem 15: The center variety of the system ( ^^ ) consists of four irreducible 
components: 

1. ViJi"^), where 4"^ = (2aio - bio,2boi - aoi), 

2. Y{4'^), where 4"^ = (aoi,6io). 
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3. V(4''^), where J^^^ = (2aoi + 601, ciio + 2&io, aoi&io - a-i2&2 -1), 

4. V{j[''^) = (/i,/2,/3,/4, /s), where /i = aj^^fes-i - a-ublo, /2 = aiofloi - 



a01^2,-l&ni- 



To our knowledge the notion of the center variety was introduced in the ht- 
erature on the center problem very recently by Zolcidek in (^91), where he also 



gives Theorem |T5|. But he writes out the fourth component in the form 

/l = /2 = 0. 

It is not precise because, using Theorem |T2|, one can see that 

v((/i,/2)) = v(ji^))uv(j('')). 

Lemma 2: The ideals I^'^^ and I^^^ are radical ideals in C[aio, aoi, c[_i2, &2,-i, ^lo, ^01] 
and C[a2o, oii, • • • , &02], respectively. 

Proof: Computing the intersection of the ideals Jk we find 

/W = Ji n J2 n J3. 



Hence /g is radical because due to Theorem |T2| Ji, J2 are prime, and according 



to Theorems |^ and ^ J3 is prime as well. 

Similarly, for quadratic system we easily check that, for 1 < i < 4, the ideals 
Jj are prime and 

4'' = nti4^^ (35) 

This yields that the ideal of focus quantities of quadratic system is a radical 
ideal. □ 

Proof: (Proposition I) V(J('')) = V{4'^^) (V(jW) = V(J^''^)) and the ideal ^'^ 
(4'^) is a radical ideal. Therefore /(^^ = 4"^ (/^ = 4'^). □ 

In (^6]), we study the following system. 

x = X- aiox"^ - aoixy - a_l3?/^ ,^^. 

y = -{y-boiy'^-bioxy-b3^^iX^). 

Using the methods of this paper, we compute the ideal of the center variety of 
(^) which turns not to be radical. Moreover, we prove the following theorem. 



Theorem 4.1: The cyclicity of the origin of the system 

ix 
is less than or equal to 5. 



ix = X — OioX^ — aQixx — a_i3X^, 
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